Stability and evolution of wave packets in strongly coupled degenerate plasmas 
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We study the nonlinear propagation of electrostatic wave packets in a collisional plasma com- 
posed of strongly coupled ions and relativistically degenerate electrons. The equilibrium of ions 
is maintained by an effective temperature associated with their strong coupling, whereas that of 
electrons is provided by the relativistic degeneracy pressure. Using a multiple scale technique, a 
(3-|-l)-dimensional coupled set of nonlinear Schrodinger-like equations with nonlocal nonlinearity 
is derived from a generalized viscoelastic hydrodynamic model. These coupled equations, which 
govern the dynamics of wave packets, are used to study the oblique modulational instability of a 
Stoke's wave train to a small plane wave perturbation. We show that the wave packets, though 
stable to the parallel modulation, becomes unstable against oblique modulations. In contrast to the 
long-wavelength carrier modes, the wave packets with short-wavelengths are shown to be stable in 
the weakly relativistic case, whereas they can be stable or unstable in the ultra-relativistic limit. 
Numerical simulation of the coupled equations reveals that a steady state solution of the wave am- 
plitude exists together with the formation of a localized structure in (2-1-1) dimensions. However, 
in the (3+l)-dimensional evolution, a Gaussian wave beam self-focuses after interaction and blows 
up in a finite time. The latter is, however, arrested when the dispersion predominates over the 
nonlinearities. This occurs when the Coulomb coupling strength is higher or a choice of obliqueness 
of modulation, or a wavelength of excitation is different. Possible application of our results to the 
interior as well as in an outer mantle of white dwarfs are discussed. 

PACS numbers: 52.27.Gr, 52.27.Aj, 52.35.Fp, 52.35. Mw 
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I. INTRODUCTION 

The nonlinear propagation of electrostatic waves in 
strongly coupled plasmas has been of considerable inter- 
est in recent years because of its possible applications 
in compact astrophysical objects (e.g. the interior of 
white dwarfs, neutron stars, the core of pre-Supernova 
stars), giant planetary interiors (e.g. Jupiter), in lab- 
oratory (e.g. ultra-cold plasmas by laser compression 
of matter) as well as in nonideal plasmas for industrial 
applications (see e.g. some recent works [ll-Hl)- White 
dwarfs have typically masses approximately between 0.07 
and 8 — lOM©, where Mq is the solar mass. It has 
been predicted that the majority of white dwarfs have 
a core of carbon-oxygen composition, which itself is sur- 
rounded by a helium layer and, for most known white 
dwarfs, by an additional hydrogen layer. However, re- 
cent research has also indicated that there may be sev- 
eral white dwarfs with large volumes (or low density) 
primarily composed of carbon with little or no hydrogen 
or helium Furthermore, Koester in a recent study 
assumed that the core may be pure C or pure O [H . The 
core of these stars is, however, extremely dense and con- 
sists of a plasma of unbounded nuclei and electrons, i.e., 
positively charged ions providing almost all the mass (in- 



' |apm isra@visva-bharati.ac.in; apmisra@gmail.com! Permanent 
address: Department of Mathematics, Siksha Bhavana, Visva- 
Bharati University, Santiniketan-731 235, India. 

f ips@tp4.rub.de; profshukla@yahoo.com| 



crtia) and the pressure, as well as electrons providing the 
pressure (restoring force) but none of the mass (inertia- 
less). In these plasma environments, (typically at high 
densities ^ 10"^^ cm~'^ and low temperatures ~ 10^ K), 
the ions form a regular lattice structure and electrons 
arc relativistically degenerate with weak electrostatic in- 
teractions d, 01- Nevertheless, in an outer mantle, plas- 
mas may be composed of non-relativistically degenerate 
electrons (at low-densities ^ 10^^ cm^"^) with strongly 
coupled ions. In both these cases, the electron Fermi 
energy becomes larger compared to both the thermal en- 
ergies ksTe^i of electrons and ions as well as the typical 
electron-ion interaction energies. Therefore, the electron 
density, to a good approximation, can be assumed to be 
uniform and unaffected by the ions. 

On the other hand, above a certain mass (M w 
IAMq), the internal pressure of white dwarfs becomes 
high enough for electrons to have sufficient momenta. In 
this case the relativistic effects become significant, and 
the degenerate equation of state P ^ r^J"^ (weakly rela- 
tivistic) that keeps the white dwarf from collapsing un- 
der its own gravityjChanges to a different form P ^ n'J^ 
(ultra-relativistic) [8[ . In the latter case, the white dwarf 
collapses under self-gravitation into a denser object such 
as a neutron star. This change of degeneracy pressure (of 
course, there is no sharp transition) will certainly change 
the spectrum of collective oscillations of electrons and 
ions, which can be used to probe the density profiles and 
other characteristics of white dwarfs. 

Recent theoretical developments have indicated that 
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strong correlations of ions significantly modify the dis- 
persion properties of collective modes as well as the char- 
acteristics of nonlinear localized structures, e.g., solitons, 
shocks in strongly coupled degenerate plasmas [HQ. A 
number of works can be found in the literature dealing 
with the linear and nonlinear properties of wave modes 
in strongly coupled dusty plasmas An addition, 

in a nonlinear regime, is the one dimensional propagation 
of wave packets and associated modulational instability 
plilin strongly coupled dusty plasmas by Veeresha et 
al 0]. Very recently, Shukla et al [l| have shown that 
the localization of ion modes is also possible in a one- 
dimensional strongly coupled relativistically degenerate 
plasma. It is, therefore, pertinent to investigate the in- 
fluence of strongly coupled ions as well as the relativistic 
degeneracy pressure of electrons on the propagation of 
multi-dimensional wave packets in coUisional electron-ion 
plasmas, such as those representative of white dwarfs. 

In contrast to (l-l-l)-dimensional evolution of envelope 
solitons, the nonlinear Schrodinger equation in multi- 
dimensions [(2-1-1) or (3-1-1)] with cubic and/ or quadratic 
nonlocal nonlinearities may no longer be integrable. In 
this case, the system often exhibits wave collapse near the 
singularity even with a wide range of initial conditions in- 
stead of stable oscillations [l2| . A collapsing wave packet 
with higher amplitudes thus self-focuses in shorter scales 
at the singularity until other physical effects intervene 
to arrest it. Nevertheless, nonlocality has been known 
to give rise to novel phenomena of generic natures. For 
example, it may promote the modulational instability in 
self-defocusing media or can arrest wave collapse in multi- 
dimensional self-focusing media. Furthermore, nonlocal 
nonlinearity may affect soliton interactions as well as may 
support the formation of stable localized structures [ij] . 
It is thus of interest to investigate the modulational in- 
stability as well as the localization of wave packets in 
strongly coupled relativistically degenerate plasmas. 



II. ANALYSIS 

A. Basic equations 

We consider the nonlinear propagation of low- 
frequency {kVn < uj < kVrc^^pi) electrostatic waves 
in an unmagnctized coUisional plasma consisting of in- 
ertial multi-charged strongly coupled classical ions and 
inertialess relativistically degenerate electrons with weak 
interparticle interactions. The dynamics of such waves 
can be described by a generalized hydrodynamic model 
which reads [ll, i, [ifl Ej 
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where n-,-, v-,- and rrij respectively denote the number 
density (with equilibrium value n^o), velocity and mass 
of different species {j = e for electrons and j ~ i for 
ions), Zi is the ion charge number, e is the elementary 
charge and Vin is the ion-neutral collision frequency. Fur- 
thermore, (j) is the electrostatic scalar potential, {Pi) 
is the electron (ion) pressure to be defined later and 
d/dt = d/dt + Vi • V is the convective derivative. In 
equilibrium, the overall charge neutrality condition gives 
neo = ZiUiQ. In Eq. ([2]), Tm is the viscoelastic relaxation 
time given by [lol | 



C + 47?/3 



1 
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where Ti is the ion temperature, ks is the Boltzmann 
constant, C + 477/3 is the coefficient of the effective ion 
viscosity in which 77 and C account respectively for the 
shear and bulk viscosity. Also, 7; is the ion adiabatic 
index and u{Ti) is a measure of the excess internal ion 
energy. Here Ti, the ion coupling parameter (to be 
presented more extensively later in the text), is a measure 
of the ratio of the Coulomb energy to the kinetic energy 
per particle. For values of F^ of the order of 1 or larger, 
correlation effects become important and the plasma is 
then called strongly coupled. The expression for u(Ti) 
can be given for k — ^ (where k, to be explained later, 
is the ratio of the interparticle distance to the ion Debye 
length) as [H [11 
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(r. < 1) , 

(1 <r, < 160), 
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(300 < r, < 2000) . 

(6) 

The relaxation time t„i represents a characteristic time 
scale to describe two classes of wave modes with fre- 
quency w, namely the hydrodynamic modes with Ci;r,„ <C 
1 and modes with ujTm 3> 1, i.e., the kinetic modes. 
The compressibility parameter ^i appearing in Eq. ([S]) 
is given by [l3| 
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9 dr, • 



(7) 



Since u{Ti) is negative for increasing values of F^, fii can 
change its sign. It has been shown that for F.^ within 
the range 1 < F^ < 10 this change of sign can cause 
the dispersion curve to turn over with the group velocity 
going to zero and then to negative values In the 

following subsections we will define few parameters that 
arc relevant in the present theory. 
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B. Degeneracy parameter 

The degeneracy parameter for a particle of species j 
can be defined as 



Xj 



^F] _ 1 /'q„2\2/3 



T, 2 



{in^f\n,X%f\ (8) 



where Xbj = h/ y/ksTjirij is the thermal do Broglie 

wavelength, Epj = kBTpj = fi^ (STr^rij)^^'^ /2mj is the 
Fermi energy and Tpj is the Fermi temperature. Thus, 
depending on the thermal energy ksTj, particles are said 
to be degenerate if the number density nj exceeds the 

quantum concentration Uqj = (2mjkBTj/h^^^ ^ /Btt^. 
Typically, for astrophysical dense plasmas, rie ^ 10^^ 
cm""^. So, Xe > 1 for Te < 10^ K. However, in metals 
(rie ~ 10^^ cm^'^) electrons are degenerate at < 10^ K. 
Thus, when the electrons form a degenerate system and 
ions are classical, Xe > 1 and Xi < 1 must be satisfied. 



C. Coupling parameter 

In order to measure the weak or strong interparticle 
interactions between electrons or ions we define two cou- 
pling parameters, namely the quantum coupling and the 
Coulomb coupling parameter. The quantum criterion of 
ideality for degenerate electrons has the form 



(9) 



where d denotes the mean interparticle distance (Wigner- 
Scitz cell radius) given by (47r/3)d'^nj = 1 and oq is 

2/3 

the Bohr radius. Since Epe rie , the parameter Fe 
decreases below one with increasing values of the elec- 
tron density (rig > 10^"^ cm~^). This implies that a de- 
generate electron plasma becomes even more ideal with 
thermal compression. Moreover, at higher densities only 
electrons represent an ideal Fermi gas whereas the ion 
component is nonidcal. Depending on the degree of its 
nonideality, one may look for ionic liquids, cellular or 
crystalline structures. Thus, one defines the criterion of 
ideality/ nonideality by the ratio between the average 
potential energy of Coulomb interaction and the mean 
thermal energy characterized by the temperature Ti for 
nondegenerate multi-charged ions as 



F,e- 



ksTid 



(10) 



where the factor k = d/Xjji measures the screening of the 
ion charge by the plasma over a distance of the ion De- 
bye length Xoi = \JkBTi/ Aima^Zfe^ . Notice here that 
in a two-component electron-ion plasma in which elec- 
trons are degenerate and ions are classical, the screening 
distance of a test charge is truly the ion Dcbye length. 



I.e. 
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where Xpe = (fcs7Ve/47rnee^) is the Thomas-Fermi 
length (screening distance by degenerate electrons) . Fur- 
thermore, plasma represents a gaseous medium for F^ < 
1, which changes to liquid-like (F^ > 10) with density 
growth or with cooling. This plasma has a short-range 
order. With further enhancement of F^, the ion sub- 
system crystallizes in the range (170 ^ F^ < 180), and for 
Ti > 180, the plasma behaves like a solid with long-range 
order. In the latter, the Coulomb force completely dom- 
inates the dynamics, and the ions arrange themselves in 
a periodic lattice structure, which minimizes the electro- 
static interaction energy. Typically, for a density ~ lO^'^ 
cm~'^ and a composition of pure O^^ in the interior of a 
white dwarf, the resulting temperature for the onset of 
crystallization is ~ 10*" K [T7| . 



D. Degenerate equation of state 

For degenerate electrons, the energy distribution is no 
longer a thermal distribution, but one governed by the 
exclusion principle, i.e., the Fermi energy which has a 
pressure associated with it. The equation of state for 
relativistically degenerate electrons is given by Q 



P. 
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3h^ 



(2s2_3) (l + s2)i/%3sinh-i(s) 
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where Pe = (Sh^rie/STT^ is the momentum of electrons 
on the Fermi surface, h {— 2'Kh) is the Planck's constant 
and s = Pe/rUeC is the nondimensional parameter. Thus, 
in the weakly relativistic (or nonrelativistic) limit (s <C 1) 
and ultra-relativistic limit (s » 1) the pressure equation 
(fT2l) reduces to two different forms: 




37r2)^/^ nl'^ = lEpene for s < 1, 



for s > 1. 



These pressures can be combined to write 



where 



Pe = K. 
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(13) 



(14) 



(15) 



(11) 



with i?5/3 = l/we and R4/3 = c. Thus, when the Fermi 
energy becomes higher than the rest energy of electrons, 
the weakly relativistic pressure equation will no longer 
be valid, and one must use the ultra-relativistic equation 
of state. The latter in the case of massive white dwarf 
stars (M w 1.4A/sun) gives rise to a lower degeneracy 
pressure than the weakly relativistic case making them 
gravitationally unstable. 



E. Pressure equation for ions 

We now consider the pressure of ions as given by [l[ 
VP, =7,M/Vn,, (16) 
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where Tif = + HiTi is the effective ion temperature 
in which appears due to electrostatic interactions be- 
tween strongly coupled ions, and is given by [isl [l9l | 



(17) 



Here, A^„„ is determined by the ion structure and corre- 
sponds to the number of nearest neighbors (e.g., in the 
crystalline state, Nnn = 12 for the fee and hep lattices, 
Nnn = 8 for the bcc lattice). Although, the parameter fii 
can be negative [cf. Eqs. © and (O] for increasing val- 
ues of Fi, may be comparable or even dominate over 
pLiTi for Ti » 1, and so the effective temperature Tif 
(> 0) is most likely due to the strong coupling of ions. 
Typically, for = 2 x 10^^ cm^^^ ^ /^)t, = W K 
and Zi = 8 (relevant for weakly relativistic regime), we 
have T, « 202, = 2 x 10^ K, /x.T, = -2.4 x 10^ K. 
Similarly, considering parameters in the ultra-relativistic 
regime, we find that is always a few orders of magni- 
tude higher than the kinetic temperature of ions, Ti. 



ne = 2 X 10^6 cm-3, Tg = AOT, = lO'^ K in the weakly 
relativistic case. In the ultra-relativistic regime, we have 
r, = 165 (Z, = 6) and = 267 {Z, = 8) for = lO^^ 
cm~^, Tf. = 2Ti = 3 X 10^ K], the hydrodynamic modes 
may not exist in the limit wr^ ^ 1 by the same reason 
as for Tm- Thus, in our coUisional plasmas we can specify 
the following two cases of interest: 

Case I: s < 1, Xe > 1, Xi < 1, < 1, > 100 
and ujTjn 3> 1. This represents the propagation of kinetic 
wave modes in a plasma, composed of weakly relativistic 
degenerate electrons with weak interactions and strongly 
coupled classical ions. In this case, the valid regimes 
for the density and temperature to be relevant, e.g., in 
an outer mantle of a white dwarf are S 10^^ cm""^, 
T,<Te< 10^ K with Z, > 1. 

Case II: s > 1, Xe > 1, < 1, < 1, > 100 and 
i-^Tm 3> 1. That is, kinetic modes may exist in a plasma 
in which degenerate electrons are ultra-relativistic with 
weak interaction and ion components form strongly cou- 
pled classical system. The parameter regimes in this case, 
to be relevant in the core of a massive white dwarf are 
He > 10^5 cm-3, T^<Te< 10* K and Z, > 1. 



F. Parameter regimes 

Our aim is to consider a fully ionized two-component 
plasma in which electrons are relativistically degenerate 
with weak interaction (almost free) and ions are multiply 
charged forming a classical system (nondegenerate) with 
strong electrostatic interaction. Since, s ^ 10~^''ne^'^, 
in the weakly relativistic limit (s <C 1) we have Ue < 
10^^ cm~'^. In this regime, the degeneracy condition for 
electrons [cf. Eq. ©] is satisfied for Tg < 10^ K. Also 
satisfied are the coupling parameters < 1. Xi ^ li 
and Ti > 1 for Zi > 1 and Ti < Tg. On the other 
hand, for ultra-relativistic degenerate electrons, we have 
Ue > 10^^ cm~^ and Tg < 10* K in order to satisfy s > 1 
and the degeneracy condition In this case, ions are 
nondegenerate (xi < 1) for Zi > 1 with a mass rrii ^ 
ZiTTiQ, where uiq = 1.66 x 10 ^^"^ g is the unit of atomic 
mass, and the conditions Xe ^ 1, Te *C 1, > 100 arc 
satisfied. 

Again, it has been shown that the viscosity coefficient 
Cvis = (C + 477/3)/riiofcsT'j in Eq. ([S]) has a wide min- 
ima (~ 1) in 1 < Fi < 10, and it tends to increase as 

oc F-/^ for Fi > 10. Also, it becomes high for F^ < 1 
Typical values of Cvis are ~ 45 for F^ = 0.1, ^ 46.4 for 
F, = 100 and ~ 118.5 for F^ = 202. Thus, (oc Cvis) 
becomes high in the weak coupling (F^ < 1) as well as in 
the high coupling (F^ > 1) regimes. So, the kinetic modes 
exist only for (F^ < 1) or (F.; ^ 1) where the condition 
^ 1 is satisfied. In the range of 1 < F^ < 10, Tm 
is typically of the order of unity, and so the kinetic con- 
dition may no longer be valid. Again, since F^ > 100 is 
more likely the case in both weakly relativistic and ultra- 
relativistic regimes [e.g. for parameters relevant to the 
conditions of white dwarfs (C^^ and O^^ compositions) 
we have F^ = 125 {Z, = 6) and F^ = 202 {Z, = 8) for 



G. Normalized system 

It is customary to normalize the physical quantities by 
redefining them in terms of new variables: Nj = nj/rijQ, 
V, = {U,V,W) = v»/Ft, T = tujp, and {X ,Y,Z) = 
(x, y, z)/Xd, where = yZ-fiksT^f /AnriioZfe'^ is the ef- 
fective Debye length, Vr — \/ jiksTif /rrii is the effective 
ion thermal speed and uipi = Vt/Xd is the ion plasma 
frequency. We then recast the basic equations ([T]), ([2]) 
and ([4]) in the following nondimensional forms. 



dT 



N,V ■ V, = 0, 



(18) 



»yv'v, + (c + ^) v(v-v,) 



(19) 



(20) 



Here d/dT = d/dT + V, • V with V = 

(77, C) = {'>l,0<^pi/niQ^ikBTif, and the constant 
is given for weakly relativistic (7 = 5/3) and ultra- 
relativistic (7 = 4/3) cases as: 



ZiTpehiTif 



for 7 = 5/3, 



pZim^c^ /"fiksTif for 7 = 4/3, 



1 /3 

where (3 — \c (STr^rieo) is the dimensionlcss param- 
eter with Ac = h/m(,c denoting the reduced Compton 
wavelength. We are then left with Eq. ([3]), which can 



(21) 
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be integrated to obtain the expression for the electron 
density N^. = n^/n^Q as 

A^e = (1 + <t>-t)^' ~ 1 + + B-ycl)^^ + a,(f>^, (22) 

where 0^ (< 1) is different for different 7, and given by 



e4>/kBTpe for 7 = 5/3, 
~ ^ ecji/PmeC^ for 7 = 4/3. 

The coefficients appearing in Eq. ([22]) are 

1 _ 2-7 

^ " ^ " 2(7-1)2' 



(23) 



(2-7)(3-27) 
6(7-1)3 ■ 



(24) 



III. 



PERTURBATION METHOD: DERIVATION 
OF COUPLED EQUATIONS 



We consider the propagation of slowly varying weakly 
nonlinear wave envelopes propagating with a group veloc- 
ity Vg oriented arbitrarily in the xy-plane to the direction 
of propagation. Then in a coordinate frame moving with 
the speed fn, the space and the time variables can be 
stretched as [H [2l| 

^ = e(.T - Vgxt), T] = e{y- Vgyt), C = ez,T = e^t, (25) 

where e is a small parameter representing the strength 
of the wave amplitude, and Vgx, Vgy are the components 
of the group velocity along the x and y axes. Since we 
are interested in the modulation of a plane wave as the 
carrier wave with the wave number k and frequency w, 
the dynamical variables can be expanded as 

00 CJO 

A = Ao + ^e" ^ 4")(C,ry,C,r)e*(''-^-*)', (26) 



ri=l 



/ — — 00 



where = 1 for Nj and = for other variables. 

Also, satisfies the reality condition A'^^i = aI"-** 

with asterisk denoting the complex conjugate. We then 
substitute the stretched variables from Eq. (f25|) and the 
expansion from Eq. (j26p into the normalized equations 
([T8 l) -p0 l) . and equate different powers of e to obtain a set 
of reduced equations as given in the following subsections. 
Until and unless mentioned we will omit the subscript 7 
in (j)^ for simplicity, but we stress that the corresponding 
results will be different for different choice of 7 = 5/3 or 
4/3. 



A. Dispersion relation 

In the lowest order of e with n = 1, Z = 1, we obtain 
the following equations for the densities and velocities in 



terms of 



a(i) 



A. 



(1) 



where A'^ = A-^ 



'^{kx,ky)cp['\ 1^W = 0, (27) 



D^k"^. Since we are considering the 



modulation of a plane wave, A^^^ are all set to zero except 
for Z = ±1. Thus, from Eq. ([27]) one readily obtains the 
following dispersion relation: 



1 



i"2 \ 2 



iujifk'^ 

1 - iUJfjn 



-0, 



(28) 



where — A^/ and fy* = C + 477/3. The dispersion 
equation (|28p has the similar form as in Ref. m for 
strongly coupled dusty plasmas. As said before, the re- 
laxation time Tm defines two characteristic time scales to 
distinguish between two classes of wave modes, namely 
'hydrodynamic modes' {toTm ^ 1) and 'kinetic modes' 
{ujTm ^ 1). However as discussed in the previous sec- 
tion, the hydrodynamic limit may not be satisfied for the 
parameter regimes in the interior of white dwarfs. So, 
considering only the kinetic limit u;r„ ^ 1, Eq. 
reduces to 



uj'^k'il 



1]- 



X2 

^P7 



(29) 



Substitution of 77* from Eq. ([5]) into Eq. further 
gives 



±k 



T, 



15 



\2 

P7 



1 



k'^^p^ 



1/2 



(30) 

For typical plasma parameters, as in Case I or Case II, the 
first term (^1) in Eq. ([30]) under the square root is larger 
than the second term (~ —0.03), i.e., the contribution 
of is always greater than that from Ti, and hence 
uj is always real for all k. However, in absence of T*, 
and since the second term may be negative, there must 
exist a critical wave number k^ below which the kinetic 
modes exist. Furthermore, Eq. ([5(1]) shows that the phase 
velocity of the carrier modes, w/fc > 1 or, a; > kVr (in 
dimensional form). Again, since the term in the square 
root is typically > 1 for ^ 1, the low- frequency (< 
ijpi) wave propagation is possible for A: < 1. Thus, in 
both the weakly and ultra-relativistic regimes as in Case 
I and II, to approaches unity for k < I. The latter also 
ensures that the wavelength is greater than the effective 
Debye length for the collective behaviors of the plasma 
not to disappear. Equation ([30]) also clears that the wave 
frequency increases with k, and it approaches the plasma 
frequency ujpi (i.e. w « 1) at fc « 0.5. For the behaviors 
of the modes we plot w vs fc as in Fig. 1. In both the 
weakly relativistic and ultra-relativistic cases, the wave 
frequency is shown to increase with increasing coupling 
parameter F^. 

We note that in the kinetic regime, the ion modes do 
not experience any viscous damping. The latter can. 
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FIG. 1. (Color online) Plot of the normalized wave frequency u vs the normalized wave number k [Eq. (|30p ] for the kinetic 
modes. The solid (blue), dashed (black) and dotted (red) lines respectively correspond to Ti = 125, 202, 512 in weakly 
relativistic case, and Vi = 165, 267, 678 in ultra-relativistic regime as indicated in the figure. The other parameter values are 
as in Case I and Case II. 



however, occur in the hydrodynamic regime. Further- 
more, the ion coupling corrections appear through u {Ti) 
whose estimates can be obtained from Eq. ([S]). Since 
u(Ti) can be negative for increasing values of F^, these 
corrections can lead to the turnover effect in which the 
frequency (and hence the group velocity) goes to zero 
and then to negative values ^M,- However, this is not 
the case for the kinetic wave modes as evident from Fig. 
1. Such effects can be strong in the 'hydrodynamic limit' 
. The term oc appears due to degeneracy pressure 
of electrons in weakly relativistic (7 = 5/3) and ultra- 
relativistic (7 = 4/3) limits. As is clear from Fig. 1 that 
the influence of this term on the dispersive curves is, how- 
ever, weaker in the ultra-relativistic limit than that in the 
weakly relativistic one. There is an additional correction 
term (cx T*) which arises due to electrostatic interaction 
of strongly coupled ions. We further note that the kinetic 
modes exist only in the strong coupling regime (F^ > 100, 
i.e. close to crystallization) where the condition a;r,„ 3> 1 
is satisfied and can be quite large. Physically, the exis- 
tence of large values of implies strong memory effects 
in the medium and hence the predominance of elastic ef- 
fects In contrast to the weakly relativistic case, the 
dispersion curves (see Fig. 1) of the ion modes for differ- 
ent values of F^ show that the wave frequency lo typically 
varies linearly with fc < 1 in the ultra-relativistic limit. 



We also find that the effect of the ion coupling parame- 
ter Vi on the kinetic modes is, however, stronger in the 
weakly relativistic limit than the ultra-relativistic case. 



B. Group velocity 

In the second order expressions for n = 2, / = 1 we 
obtain the corrections for iVg\ , A^j\ etc. in terms of 
(p\ and cpl . After eliminating those iVe\ ' ^ii etc., we 
obtain a resulting equation in which the coefficient of 
vanishes by the dispersion relation, and the coefficients 
of dcj}^^^ /d£^ and dfj)^^^ /drj are set to zero to obtain the 
group velocity components (compatibility condition) as 

duj 2k 

(31) 

where 

^ = 2u:+—^-—^. (32) 

(1 - iujT„-,y 

The expression for the group velocity can further 
be simplified in the limit of ujt„i 3> 1 and us- 
ing Eq. (O. For wr™ ^ 1, we have ui = 
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FIG. 2. (Color online) The normalized group velocity dispersion Vg [Eq. (|33[l] is plotted against the normalized wave number 
k for the kinetic modes. The solid (blue), dashed (black) and dotted (red) lines are for different values of F,; as in Fig. 1. 



(2/sVw) [u}^/P - ii-q* /fm)/2ujfm] ■ Since ?7*/f„ < 1 for 
Ti ^ 1; also uJTm 3> 1 (already assumed) and oj/k > 1 
from the dispersion relation, the imaginary part of Cj can 
be neglected compared to the real part. Thus, we obtain 



k 



T^f 



— — 1 H u + -^-^ 



(33) 



From Fig. 2 we find that in the weakly relativistic 
limit [Fig. 2(a)], the modes with short-wavelengths have 
greater group velocity than those with long-wavelengths, 
and this group velocity is always smaller than the effec- 
tive ion thermal speed as well as the phase speed of the 
carrier wave. Also, in this case Vg increases as the wave 
number k (< 1) increases. However, as k approaches 
(i.e. in the long wavelength limit), Vg tends to a steady 
state value. This anomalous group velocity dispersion is 
an important condition for the modulation instability of 
wave packets. In contrast to the weakly relativistic case, 
the wave modes with long-wavelengths propagate with 
higher group velocity than those with short-wavelengths 
in the ultra-relativistic limit [Fig. 2(b)]. Also, the group 
velocity in this case is greater than the ion thermal speed 

I 



and the phase speed except for some higher values of 
Ti ^ IC*. The latter may not be relevant in in the inte- 
rior of white dwarfs as values of Zi are not so large there. 
Furthermore, in the ultra-relativistic limit, the group ve- 
locity tends to decrease with increasing values of < 1, 
and there exists a critical value kc < 1 of k above which 
the behaviors of Vg remain almost the same even with dif- 
ferent values of the ion coupling parameter Ti. In both 
these weakly relativistic and ultra-relativistic cases, the 
higher the coupling parameter Ti, the lower is the group 
velocity of dispersion. 



C. Coupled equations 

Considering the zeroth harmonic modes for n = 2,3, 
I — 0, which appear due to the nonlinear self-interaction 
of the carrier waves, we obtain corrections for the den- 
sities and velocities (see for details, Appendix-B). These 
expressions are then used to eliminate the variables to 
obtain the following equation for 



g2 q2 q2 q2 \ 



g2 q2 q2 q2 \ ^ 



(34) 
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where the coefRcients are given by 



Ujkr 



(35) 
(36) 



5*3 = 2B^ 



Sa 



(37) 



r 



We ought to mention that in the coefficient of for n = 
2, / = 0, one usually obtains from the momentum balance 
equation for ions the relation S'inl'/'i^^P = 0. Since 
is of the order of , should be at least of the order of 
e, and so it will contribute to the coefficient of e'^ of the 
momentum equation for r? = 3, / = 1. 

Next, the second order harmonic modes for n = 2, Z = 2 
are obtained as 



U. 



(2) 



A 



22 [ 



(2) 



(2) 



e2 



A. 



7'^2 



B^ 



(39) 



(2) 
i2 



(2) 
e2 



a(2) 



7V^2 



4^) 



D 



22 1 



(40) 



where the coefficients arc given in Appendix C. Finally, 
for n = 3,^ = 1 we obtain expressions for the third- 



order first harmonic modes. Here the coefficients of 



and 



a(2) 



vanish by the dispersion relation. Thus, we ob- 



tain after a few steps the following nonlocal nonlinear 
I 



Schrodinger equation: 



dr 



Pi 



92 



6W 



d^dr] 
(41) 



In deriving Eq. (|4ip , we have neglected a small contri- 
bution from the term {kxUlf^ '^1^'' compared 

to dcjy'i'' /dr. Because, from the dispersion relation we 
find that the phase speed of the wave is larger than the 



(38) '^'^^ thermal speed. Also, the modes Uq"^^ and Vq^' , which 
arise due to the mean motion of ions, can even be smaller 

than the ion thermal speed, since they appear as small 

(2) (2) 

(higher order of e) corrections. So, w ^ k^U^ ' + kyV^ ' 
can be a good approximation to the present case. 

The terms with coefficients Pi , P2 appear due to wave 
group dispersion, those with P3 and P4 are respectively 
due to the departure from the two-dimensional formalism 
and due to the oblique modulation. The cubic nonlinear- 
ity (Kerr) is due to the carrier wave self-interaction origi- 
nating from the zeroth harmonic modes (or slow modes) , 
and the nonlocal nonlinear (quadratic) one appears due 
to the coupling between the dynamical field associated 
with the first harmonic (with a 'cascade effect' from the 
second harmonic) and a static field generated due to the 
mean motion (zeroth harmonic) in plasmas. 

Now, Eq. (|4T|) is coupled to Eq. (p4|) . and we recast 



r(2) 



a(i) 



the two equations as [considering now that 
dependent on 7, i.e. rewriting $^ 
for weakly relativistic (7 
(7 = 4/3) cases] follows: 



and 



are 



5/3) and ultra-relativistic 



p5'*7 



P2^ 



3 Q(^2 



(42) 



R2^ 



^3 ^TT^ r it4 - 



d^drj 



^3^ 



d£,dr] ' 



(43) 



in which the coefficients are different for different values are weakly relativistic or ultra-relativistic) and are given 
of 7 = 5/3,4/3 (according as the degenerate electrons 



9 



in Appendix D. 



IV. MODULATIONAL INSTABILITY 

The modulation of slowly varying wave amplitudes 
may occur due to, e.g., parametric wave coupling, nonlin- 
ear interaction between high- and low-frequency modes 
or, self-interaction of the carrier wave modes. These phe- 
nomena, however, have relevance with the MI, which con- 
stitutes one of the most fundamental effects associated 
with the wave propagation in nonlinear media. Such in- 
stability basically signifies the exponential growth or de- 
cay of a small perturbation of the wave as it propagates 
in plasmas. The gain leads to amplification of sidebands, 
which break up the otherwise uniform wave and lead to 
energy localization via the formation of localized struc- 
tures. Thus, the MI may act as a precursor for the for- 
mation of bright envelope solitons. On the other hand, 
the formation of dark solitons requires the absence of MI 
in the constant intensity background. The MI can, how- 
ever, be affected by the obliqueness of modulation, the 
electron degeneracy as well as the strong coupling effects 
of ions which wc will discuss shortly. 

Let us consider the modulation of a plane wave solu- 
tion of Eqs. (|42|) and (|43|) of the form (omitting the sub- 
script 7 once again, for simplicity) $ = $oe~'^°^: ^' = 
so that = — (3i$g — iQs, where $o is a constant. 
Here the choice of ^ is immaterial as the stability condi- 
tion does not depend on it, and the solution for $ is not 
unique. We then modulate the wave amplitude as a plane 
wave perturbation with frequency n and wave number 

K, i.e., $ = ($0 + $ie*K R-iOr _^ ^^g-iK R-hiOr-J g-iOor 

and * = ^^e*^-^-'^^^ + xi,^e-'^-^+'^'' , where $1,2, 
^'1^2 are all real constants, K = {Ki, K2, K3) and R = 
(^,77,0. Looking for the nonzero solution of the small 
amplitude perturbations, we obtain from Eqs. (|42|) and 
P5)) the following dispersion relation for the modulated 
wave packet 



ft 



2$2 



/i 



^2/3 



(44) 



where fi , /2 and /a are given by 

/l,2,3 = iPl,RuSi)Kf + {P2,R2,S2)Kl 
-f (P3, i?3, 53)At + (^4, i?4, Si)KiK2. 



(45) 



In the limit of uJTm ^ 1, all the coefficients in Eqs. 
and (|43|) are real and so, Eq. (|44)) can be rewritten as 



(46) 



where is the critical wave number given by 

2(J)2 , s 

Kl = ^(Qif2 + Q2h){l + al + al). (47) 
/1/2 ^ ^ 



Here, /i,2,3 = fia-.^l^i ^nd ai,2 = A'2,3/^1 define the 
obliqueness of perturbation of the wave numbers. As 
we will see later, those perturbations can change the 
stable and unstable regions for wave packets. Equa- 
tion p6)) shows that the MI sets in for K < Kc, and 
the right-hand side of Eq. (|47| is positive, i.e. T = 

{QiIi + Q2./3^ /./1/2 > 0. In this case the perturba- 
tions grow or decay exponentially during propagation of 
waves. On the other hand, for K > A'c, the wave packet 
is said to be stable under the modulation. The instability 
growth or decay rate (letting VL = iFinst) can be obtained 
as 



1 



1, 



where the maximum value is achieved at K 
and is given by 



-,(max) 



-f (Q1./2 + O2/; 



(48) 
Kc/V2, 



(49) 



We numerically investigate the stable and unstable re- 
gions for the plane wave solution relying on the above 
condition and with different coupling parameters 
as applicable for weakly relativistic as well as ultra- 
relativistic regimes (c/. Case I and Case II). We consider 
the parameter values which are relevant to the conditions 
of white dwarfs. For example, in the weakly relativistic 
case. Tie = 2 X lO^^ cm^^^ = iOT, = lO"^ K and vary 
Z,=6 (Fj = 125), 8 (F, = 202). In the uhra-relativistic 
case, we take = lO^^ cm^^, Tg = 2T, = 3 x 10^ K 
with Z, = 6 (F,; = 165) and 8 (F, = 267) for C^^ and O^^ 
compositions. Due to highly efficient thermal conductiv- 
ity properties of the degenerate electron gas, the interior 
of a white dwarf is nearly isothermal with a temperature 
~ 10^ or 10® K. 

Since K represents the wave number of perturbation 
it should not be higher than the carrier wave number k. 
However, ai and 02 can even be larger than unity, and 
thereby can change the sign of T, and hence the stable 
and unstable regions. For example, for a fixed F^ = 125 
as in Fig. 3(a), the wave is stable in the regimes, namely 
< ai < 0.5 and 02 > 0; < ai < 0.6 and 02 > 0.2; 
< ai < 0.7 and a2 > 0.4 etc. In the other regimes of ai 
and a2, the wave can be unstable. Fig. 3 shows different 
stable (white or blank) and unstable (colored or shaded) 
regions in which panels (a), (b) are for the weakly rela- 
tivistic case with different F^ = 125 and 202 respectively, 
and panels (c), (d) for F^ = 165 and 267 in the ultra- 
relativistic regime. From Figs. 3(a) and (b) we find 
that when the ion coupling strength is relatively small 
(e.g. Fi = 125 in the liquid state), the long wavelength 
kinetic modes exhibit instability against an oblique mod- 
ulational perturbation in strongly coupled plasmas with 
weakly relativistic degenerate electrons (e.g. in an outer 
mantle of white dwarfs). However, as F^ increases (e.g. 
Fi = 202 in the crystallized state), the instability do- 
main for the wave number expands, and the system can 
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FIG. 3. (Color online) The stable (blank or white) and unstable (shaded or colored) regions are shown in the plane of 
normalized wave number k and the obliqueness of modulation a, corresponding to the conditions T < and T > respectively 
as in the text. Subplots (a), (b) are for Fi = 125, 202 (weakly reiativistic), and (c), (d) (Ultra-relativistic) are for F; — 165, 
267 respectively. Other parameter values are as in Case I and Case II, and K = 0.001, ai = 2, 02 = 0.1, 0o ~ 0.01. 



exhibit oblique modulational instability for a wide range 
of values of k. From Figs. 3(a) and (b) we can also con- 
clude that the plane waves during propagation in plas- 
mas, e.g., in an outer mantle of white dwarf, exhibit sta- 
bility against small perturbations and for a wide range of 
values of a (where a is the angle, the wave vector k makes 
with the X-axis) as well as the wavelength. Furthermore, 
when the modulation takes place along the direction of 
propagation of waves (i.e. for a ~ 0) in strongly cou- 
pled plasmas with weakly or ultra-relativistic degenerate 
electrons, the system always shows stability regardless of 
the values of F^. From Figs. 3(c) and (d) we find that 
the instability of plane waves in strongly coupled plas- 
mas with ultra-relativistic degenerate electrons (e.g. in 
the interior of white dwarfs), depends strongly on the 
scale length of excitation of the wave modes, the angle 
of modulation, the ion coupling parameter as well as the 
degeneracy pressure of electrons. Thus, strong coupling 
of ions (Ti > 100) as well as the ultra-relativistic degen- 
eracy pressure of electrons favor the instability of plane 
waves against the oblique modulation. 

The instability rate Finst is calculated and the devel- 
opments of the instability growth or decay are shown in 
Fig. 4. We find that the maximum value of the insta- 
bility rate typically depends on the cubic as well as the 
nonlocal nonlinearity through the coefficients Qi and Q2 



respectively. It is clear that the electron degeneracy has 
an important effect for the change of sign of the insta- 
bility rate Finst- For the waves propagating in strongly 
coupled plasmas with weakly reiativistic degenerate elec- 
trons, the plane wave perturbation grows exponentially 
with the wave number leading to the instability growth 
as in Fig. 4 (upper panel). However, a lower value of 
Fi (when ions are in the liquid state) tends to suppress 
the instability, decreasing the growth rate at lower value 
of the wave number k. This means that the plane waves 
in plasmas with crystallized ions and weakly reiativis- 
tic degenerate electrons exhibit long wave MI. On the 
other hand, for the propagation of waves in strongly cou- 
pled plasmas with ultra-relativistic degenerate electrons, 
the instability rate becomes negative by the electron de- 
generacy. In this case, the contribution from the non- 
local nonlinearity becomes larger than the cubic nonlin- 
ear term leading to exponential decay of the perturba- 
tions, and the decay rate can not be suppressed by the 
strong coupling of ions, i.e., whenever ions are in liquid 
or crystallized states. Thus, while the growth rate can 
be suppressed (i.e. Finst < 1) by lowering the coupling 
parameter F^ in the range 100 < F^ < 200 (in the weakly 
reiativistic case, see the upper panel of Fig. 4) with a cut- 
off at K < 1, the decay rate remains unbounded within 
K < 1, and can not be controlled with a considerable 
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FIG. 4. (Color online) The growth (case of weakly relativistic, see upper panel) and decay (case of ultra-relativistic, see lower 
panel) rates [Eq. (|48|) ] are shown at a = 2 and a = 1.7 respectively. The solid (blue), dashed (black) and dotted (red) lines 
are for different values of Ti as in Fig. 1. The other parameters are as in Case I and Case II, and k = 0.2, qi — 2, a2 ~ 0.1, 
(j>o = 0.01. 



range of values of Fj. As an example, the dotted line 
in the lower panel of Fig. 4 has the cut-off at K = 7.4 
for Ti — 512 corresponding to Zi = 14. The other lines 
(solid and dashed) corresponding to F^ = 125, 202 have 
cut-offs at higher K > 7.4. 



V. EVOLUTION OF WAVE PACKETS 

In order to examine a long-time evolution of the wave 
packets, we perform a numerical simulation of Eqs. (|42p 
and (031) in the weakly relativistic regimes. To this 
end, we discretize the derivatives using a finite difference 
scheme, and use the domain size 25 < ^, < 25 with 
150 X 150 grid points and time step dr = 10~^ for (2-1-1) 
dimensional evolution. In the case of (3-|-l)-dimensional 
evolution, we use the domain size 15 < ^,77, C < 15 with 
60 X 60 X 60 grid points and time step dr = 10~^. A sym- 
metric Gaussian wave beam is chosen as an initial condi- 
tion of the form ~ ^2//7ra6exp(— ^^/a^ — if /b"^) in 
(2+1) dimensions [similar for (3+1) dimensions], where / 
is the wave action. Notwithstanding, we have made sim- 
ulations with other initial beam profiles, however, the 
qualitative results remain almost the same. Similar anal- 
ysis can also be done in the case of ultra-relativistic limit 
(7 = 4/3), however, one has to be careful about the coef- 



ficients which become larger in magnitude, and one might 
have to rescalc those in order to perform numerical anal- 
ysis with smaller step size for the space and/or time. 



A. (2+l)-dimensional evolution 

We neglect the z-dependence of the physical quantities. 
Then the coefficients P3, -R3 and 5*3 are all zero, and Eqs. 
(|35|) and P5)) reduce to Davcy-Stewartson-like equations 
in a more generalized form. The latter, in particular, has 
been inves tig ated by a number of authors (see, e.g. some 
Refs. in [l3 . [2l|) not only in the context of plasmas, 
but also in some other nonlinear media. Such equations 
may eventually give rise to unstable solutions other than 
localization, or the formation of a singularity at which 
the solution blows-up with higher amplitudes in a shorter 
scale. This means that system's validity breaks down 
near the singularity implying that an additional physical 
mechanism might be necessary in order to arrest such 
blow-up. Here we show that no singularity is formed for 
a wide range of parameters appropriate for the model, 
i.e., the model that we have considered is self-consistent 
in (2+1) dimensions, and does not give any unbounded 
solution. 

Figure 5 shows the time evolution of the wave ampli- 
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FIG. 5. (Color online) Time evolution of [numerical solution of Eqs. (|42p and (|43[l ] the wave amplitude (7 = 5/3) for 
different Ti. The solid (red) and dashed (blue) (see the inset) lines represent stable oscillations, and correspond to = 0.1 
and Vin = respectively (with k = 0.5, a = 2, Fi = 125). The dotted line shows irregular oscillations for a different set of 
parameters, i.e., Vin = 0, k — 0.2, a = 1.6 and Fi — 202. 



tudc with two difFcrcnt values of the coupling parameter 
Ti = 125 (blue or dashed line) and 202 (black dotted 
line) with (solid or red line) and without (dotted and 
dashed lines) collisional effects. The solid or red line rep- 
resents the curve with weak collisional effects and with 
the same = 125 as the dashed line. A steady state 
oscillation (see the inset), (even when the wave action 
is well above the critical power) is found, to be called 
driven damped oscillations. Usually, due to the the col- 
lisional term, which acts like a frictional force, the wave 
amplitude dies away. However, in this case, the nonlocal 
nonlinearity always feeds the energy into the system so as 
to offset the frictional losses. Since the collision frequency 
is smaller than the ion plasma frequency and Um ~ e, i.e., 
Q3 ~ e/2, higher values of Din are inadmissible, and we 
still have a stable solution. 

For a different set of values of the coefficients Pi, P2, 
etc. with a higher F^ — 202, the wave amplitude oscil- 
lates in an irregular manner (see the dotted line in Fig. 
5), implying that the perturbations are not stable, and 
the wave packets may not be localized. Because, in con- 
trast to the case of the solid or dashed curves (in which 
Pi, P2, P4 ~ Qi and Q2 < Qi), the dispersion coeffi- 
cients become larger than the nonlinear terms, and either 
the cubic or quadratic nonlinearity is not enough to bal- 
ance the higher dispersive effects. Thus, the localization 



of wave packets may be possible in plasmas with ions 
in liquid state and degenerate electrons are weakly rela- 
tivistic. However, for plasmas in which ions are in crys- 
tallized state and degenerate electrons weakly or ultra- 
relativistic, the wave packets may not be localized with 
higher values of F^ ~ 200. This is expected as we have 
seen in the previous section that the MI growth rate of 
plane waves can be suppressed by lowering the coupling 
parameter F^ instead of its higher values in the weakly 
rclativistic strongly coupled regimes. The higher values 
of Ti gives rise to the enhancement of the growth rate 
leading to an exponential growth of perturbation with 
no cut-offs at a finite fc < 1 . 

Figure 6 shows examples of time evolution of wave 
packets in the case of stable oscillations of wave ampli- 
tude. We find that the wave is localized and propagates 
with a permanent profile. The parameter values are con- 
sidered as the same as for the solid line in Fig. 5. Initially, 
the wave amplitude decays (and the symmetric Gaussian 
form breaks down) into different shape until the wave 
gets modulated. The wave amplitude then starts grow- 
ing until the maximum modulation is achieved. After 
some time, the wave gets stabilized and propagates with 
a permanent profile due to nice balance of the dispersion 
and nonlinearity. Thus, we conclude that the localiza- 
tion of wave packets having kinetic ion modes as carrier 
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FIG. 6. (Color online) Contour plots of a numerical solution of Eqs. (|42p and (|43[) at different times as in tfie figure. The 
initial wave form decays due to dispersion until it gets modulated. The amplitude then grows until it takes the maximum 
modulation, and then propagates with a permanent profile. The parameter values are the same as for the stable solution, i.e., 
the solid line in Fig. 5. 



modes in plasmas where strongly coupled ions arc in liq- 
uid state and degenerate electrons are weakly rclativistic 
such as those in an outer mantle of white dwarfs, is pos- 
sible through the MI. These packets propagate in such 
plasmas with a permanent profile for a long time. On 
the other hand, plasmas with crystallized ions and rela- 
tivistically degenerate electrons, such as those in the core 
of white dwarfs, can not support the propagation of such 
wave packets with a permanent profile. 



B. (3-|-l)-dimensional evolution 

We examine the effects of additional dispersive and 
nonlinear effects that are due to the z-dependence of the 
physical variables. We find that a singularity is formed 
at which the wave amplitude blows up in a short time 
(see the dash-dotted line in Fig. 7) for the same set 
of parameters as the case of a steady state solution in 
(2-|-l)-dimension [see dashed (blue) line in the inset of 
Fig. 5]. In this case, the coefficients are Pi = 1.59, 
P2 = 1.3, F3 = 1.68, P4 = 0.9, Qi = -1.4, Q2 = 0.35, 
Qs = 0.0, Ri = -1.8, R2 = -1.16, i?3 = -1.97, Ri = 
-0.74, Si = 0.75, 5*2 = -2.1, S3 = 1.5 and S4 = -8.0. 
That is, additional effects due to the dispersion (P3, R3) 
and nonlinearity {S3) are responsible for the formation of 



singularity. It turns out that the validity of Eqs. and 
(|43p in (3+1) dimensions breaks down near this singular 
point. Since a physical quantity can not be infinite, some 
of these dispersive and/or nonlinear effects, which were 
initially small, become important in order to prevent such 
collapse. 

We find that when the dispersive effects are more pro- 
nounced than the nonlinearities, which occurs for a differ- 
ent set of parameters, i.e. of k and/or the obliqueness a 
or by increasing the coupling parameter [e.g. k = 0.2, 
a = 2, Ti = 125 (dashed line in Fig. 7) or fc = 0.2, 
a = 1.6, Ti = 202 (solid line in Fig. 7)], the collapse 
is arrested with an unstable oscillation. We do not find 
any steady state solution for a wide range of parameters. 
The latter, however, requires further investigation for a 
conclusive evidence, and is limited to the present study. 
The weak ion-neutral collision has no effect other than a 
weak damping of oscillations (not shown in the figure) of 
the wave amplitude. 



VI. DISCUSSION AND CONCLUSION 

Results of the previous sections demonstrate that plas- 
mas with strongly coupled ions and weakly or ultra- 
relativistically degenerate electrons can support the exci- 
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FIG. 7. (Color online) Time evolution of [numerical solution of Eqs. (|42p and (|43p ] the wave amplitude $7 (7 = 5/3) for 
i^in ~ 0. The dash-dotted line shows that the solution blows-up in a finite time corresponding to the same parameters, i.e., 
k = 0.5, Q = 2 and Ti = 125 as for the dashed line (stable solution) in Fig. 5. The solid (blue) and dashed (red) lines show that 
the collapse is arrested with unstable or irregular oscillations (bounded) corresponding to fc = 0.2, a = 1.6, — 202 (same 
parameters as for the dotted line in Fig. 5), and a different set, i.e., k = 0.2, a = 2, Ti — 125 respectively. The inset shows a 
profile for vp^, corresponding to the dashed (red) line for No stable solution is found. 



tation of low-frequency kinetic ion wave modes. The lat- 
ter with short-wavelengths travel with a higher or lower 
group velocity than those with longer wavelengths de- 
pending on whether the electrons are weakly relativis- 
tic or ultra-rclativistically degenerate. This anomalous 
group velocity dispersion is one of the most important 
conditions for the modulational instability of wave pack- 
ets in nonlinear media. We have discussed the parameter 
regimes in two cases in which the present model is valid. 
These regimes are, in particular, representative of the 
interior (e.g., carbon-oxygen composition) or an outer 
mantle of white dwarfs. 

In the present model, the equilibrium of electrons has 
been considered to be maintained by two pressure equa- 
tions as per Chandrasekhar Q , whereas that of ions are 
associated with the strong coupling effects [l^, [l^ . We 
find that the dominant contribution of the effective ion 
temperature is mainly due to the strong interaction of 
ions. In absence of the latter, the wave modes propagate 
with wavelengths above a critical value, and its contribu- 
tion in the wave group dispersion is higher than the ki- 
netic temperature of ions. Since the relaxation time has 
been found to increase highly with the ion coupling pa- 
rameter Ti > 10 [2^, only kinetic wave modes {ujTm ^ 1) 



exist in plasmas where ions form a liquid state or a crys- 
tallized one. In the latter, the elasticity will dominate 
over the viscosity effects. The hydrodynamic modes, on 
the other hand, may be relevant in other regimes where 
Ti is not so large, e.g., in 1 < < 100. In this case, 
the medium behaves like a liquid where the viscosity ef- 
fects become important, and the wave frequency to has to 
be well below the ion plasma frequency in order to sat- 
isfy the hydrodynamic limit. Also, for lower values of uj 
where the dispersion is weak, the soliton (e.g. Korteweg- 
de Vries soliton) formation of the carrier waves is a lower 
order process than the modulational instability of the 
wave envelopes discussed here. 

We have studied the modulation of a waveform with 
constant power by perturbing its amplitude with a plane 
wave disturbance in multi-dimensional form. It shows 
that the obliqueness of modulation in the xy plane desta- 
bilizes the wave packet whereas it remains stable un- 
der the parallel modulation {a = 0) in both the weakly 
relativistic and ultra-relativistic regimes. Furthermore, 
in contrast to (l-l-l)-dimension, the wave number of 
perturbations in multi-dimensional propagation, changes 
the domain of stable and unstable wave numbers sig- 
nificantly. This implies that what is stable in (1-f-l)- 
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dimension could become unstable in a more general situ- 
ation. For the parameters as in Case I, the wave packet is 
shown to be stable within a wide range of the carrier wave 
number k and the obliqueness parameter a. The region of 
instability may be increased in the ka plane with higher 
values of Ti. In these higher values, the growth rate of in- 
stability remains positive and high, which, however, can 
be suppressed by a lower value of Ti with a lower cut-off 
of the wave number of perturbation K < 1. On the other 
hand, the results corresponding to Case II show that the 
wave envelope is modulational unstable for a wide range 
of k and a. In this case, though the decay rate can be 
lowered by increasing F,;, however, it remains outside the 
considerable range of K < I. Thus, for the propagation 
of wave packets in plasmas with ions in liquid or crystal- 
lized states as in the core of white dwarfs, the decay rate 
of instability may not be controlled. However, in plasmas 
with liquid state of ions and non-relativistic degenerate 
electrons (as in an outer mantle of white dwarfs), the 
waves are mostly stable to modulational perturbation, 
and the instability can even be suppressed with a lower 
value of Ti. This implies that the localization of wave 
packets having kinetic ion modes as the carrier waves as- 
sociated with the MI may be possible in plasmas (e.g. in 
an outer mantle of white dwarfs) in which ions are in liq- 
uid state and degenerate electrons are weakly relativistic. 
Since the growth or decay rate of instability can not be 
suppressed by higher values of F.^, the MI does not give 
rise to localization of wave envelopes in plasmas where 
ions arc in crystallized states and electrons arc relativis- 
tically degenerate (e.g. in the core of white dwarfs). 

Next, we examine whether the localization of wave 
packets is possible though the MI in some regimes where 
the instability growth rate is suppressed, and in some 
other regimes in which the time evolution of the wave 
amplitude exhibit irregular oscillations leading to delo- 
calization of the wave envelopes. The time evolution of 
the (2-f l)-dimensional wave packets shows that the wave 
amplitude stabilizes for a long time with a lower value of 
the coupling parameter F^ ^ 125. However, for a higher 
value of Fi ^ 202, it shows irregular oscillations (see Fig. 
5), i.e., unstable. In the case of stable wave oscillations, 
we find that for an initial Gaussian wave beam, the am- 
plitude initially decays and the beam looses its shape be- 
fore it gets modulated by the nonlinear self-interactions. 
The amplitude then starts increasing until the maximum 
modulation is reached. As the time progresses, the wave 
amplitude reaches a steady state value, and the initial 
beam transforms into another with a permanent profile 
(see Fig. 6). We have included a phenomenological ion- 
neutral collision term which appears in the third-order 
corrections as a linear term with the first order pertur- 
bation, and remains weaker than the nonlinear (Kerr and 
nonlocal) contributions. It has no effect on the stability 
or instability of wave modes under modulation, rather it 
changes the oscillation pattern with a higher amplitude 
like a steady state damped harmonic oscillation (see solid 
line in Fig. 5). The frictional force due to the collision. 



which usually damps the wave amplitude, does not make 
so as the nonlocal nonlinearity feeds up the sufficient en- 
ergy in order to offset the decay. 

On the other hand, the results in the (3+1)- 
dimensional evolution indicate that a wave singularity 
can be formed leading to wave collapse in a finite time 
when the dispersive effects are not in balance with the 
nonlinearities, or the additional dispersion (due to z co- 
ordinate) dominates over the group velocity dispersion. 
This collapse can, however, be arrested for a different 
choice of k and/or a, or by increasing the coupling 
strength F^ where the dispersion is more pronounced than 
the nonlinearities. The numerical simulation reveals that 
the present model in (3-1-1) dimensions does not support 
the formation of localized structure with stable oscilla- 
tions. In order to obtain such steady state solution, one 
way could be to consider Vg in an arbitrary direction of 
space or along a fixed axis. However, one needs further 
investigation in order to confirm it, which is limited to 
the present study. 

To conclude, in strongly coupled plasmas with rel- 
ativistically degenerate electrons, one should carefully 
consider the parameter regimes as discussed in Case I 
and Case II. In high coupling regimes (F^ ^ 1), the effec- 
tive temperature associated with the strong electrostatic 
interactions of ions is much more pronounced than the ki- 
netic (Ti) one. The modulation of plane waves and their 
nonlinear evolution depend strongly on the system pa- 
rameters including the obliqueness a, the Coulomb cou- 
pling Ti as well as the extra dimension (geometry) of 
the system. We find that the localization of wave pack- 
ets in (2-1-1 )-dimcnsions having kinetic ion modes as car- 
rier waves is possible in plasmas with ions in liquid state 
and degenerate electrons are weakly relativistic. In other 
plasma regimes, e.g., in the core of white dwarfs where 
ions arc in crystallized state and degenerate electrons are 
ultra-rclativistic, the MI of plane waves gives rise an un- 
bounded growth or decay of instability leading to delocal- 
ization of wave envelopes. Since Ti can be controlled by 
heating or cooling the ion components, it would be inter- 
esting to look for the excitation of ion wave modes, and 
their localization as wave packets in laboratory experi- 
ments. Our model is more general than those available 
in the literature [l|, Q , and could be applicable to other 
nonidcal systems, e.g. metal plasmas. 
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APPENDIX A: SECOND ORDER FIRST 
HARMONIC MODES 



For n = 2, / = 1 we have 



a(1)' 



(50) 



= kU['^ - ^2D. I + k,,^ I , (51) 



1 — lUJTr, 



(54) 



— Ci2[x,y)^l2{y,x) — C'l2(y -^12 , (55) 



(auf, "ur,) = Bi2{Ei2y, Di2y) - Ai2y{Di2x,Ei2x), (56) 



{av^,avrj) — Bi2{Dl2x,Ei2x) — Ai2x{Ei2y,Di2y), (57) 



a 



where 

a = i?^2 ^ A\2xAi2y^ Bi2 



(52) 



2 . . o _ (C + ^/3)fc.fcy^ ^gg^ 



(^12 — — 



Cl2{x,y) = (£'7 + -^7) ' (58) 



w + irjk^ /{\ — iujTm) 



1 — iujTr, 



(59) 



Dl2{x,y) = TJf. — ^ L(l - i2a;f„)Wg(^,j,) + i2T]* k(^x^y)\ - [1 - if„, (w + k(^x.,y)Vgix^y))\ 



+1 



{C + fj/i)K^^ukl^ 



^) on i.2 



(60) 



ujK^k 



{x,y) 



P(l - iLOf^) 

.(C + ^/3) A'2^fc^fc.y 



fc2(l - iujTm) 



[(1 - z2u;r„)ug(j^,2,) + i2ryfc(j^^^)] + i- 



+ 7/ '^^'1' {x,y)'^g{y,x) 

1 - iwf,„ 



(61) 



r 



APPENDIX B: ZEROTH HARMONIC MODES 



For / = we obtain the zeroth harmonic modes from 
the coefficients of and C' in terms of ^g^' and 



AWq^^^ {A^ + 



(2) 



^eo^ = = + 2B,\^Y'{', (62) 



(1)|2 



7) Q(- -7-7 ' 

where A = Vgxd^^ + "^gyd-q and iJ^.y, C'^.y are given by 



(66) 



(2) m/(2) 



(2) 



d-q dC, 



Bx.y = 2B., + D^Kl ' - • • 



■a::!; I ^g{x.y) 



/c2 



AyJ'^ = (A^ + i?7)^^ + Cy "''?. ' + B; 



/c2 



Af/o^'^ = (A^ + i?^) 



(63) 



J2) 



^0 .^Mi^+c^''^^ 



(1)|2 



a(2) 



(64) 

a(1)|2 



_ 0^ 



Cx.v — ^ n kx,y 1 ^g{y,x) 



ujk„ 



fc2 



(67) 



(68) 



5102^ P P ^'^'^ operator A once in Eq. and ehminate 



5r7 



a^ ^ 9?7 



Uo^\ Vo^\ ^0^^ by using Eqs. 1^ and dM!)-® 



(65) to obtain Eq. 
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APPENDIX C: COEFFICIENTS OF SECOND 
ORDER SECOND HARMONIC MODES 



Ai (C + f]/3) k^ky 
1 — i2ujfm 



(70) 



The coefRcients of the second order harmonic modes 
for n = 2, Z = 2 are given as fohows: 



Jx,y 



= UJ + 



2i [fyP + (C + ^/3) fcgj 
1 — i2u)f„, 



h,2 = -K^kx,y 



(69) 



lUJTr, 



1 — i2uJTrr 



fc2 



(71) 



A22 = ^ [2kxB^ - 01 - B220xy + 2k,D22 (A + ^7(1 + 4fc'))] , B22 = ~ ^2), 



(72) 



D22 = ^ [iok^iB^ - K^){Oly - 0.0y) - {Oi - 2k,B^){2kl0y - k^kyOxy) + (92 - 2kyB^)ik%y - 2k,kye,)] , (73) 
^ = {A^+ 4k^D^) [AklikyOxy - kx9y) - wA:,(Cy - ^OxOy) - 2k^k%] - AD^k^iklOy + kp^ - k,ky9xy), (74) 



i?i = 2 [{A-, + Ak^D^){kl - LoOx) + D-,kl] [fc,(02 - 2B^ky) + 2uj9,y{B^ - K^)] , 



(75) 



7?2 - [(A + Ak^D^){2k.,ky - c^0„,) + 2D^k,ky] [k^iOi - 2S^fc,) + 2Lu0,iB^ - K^)] , (76) 

I 



APPENDIX D: COEFFICIENTS OF EQ. (|42|) 
Here we give the coefficients Pi, P2 etc. as follows: 



k^ 

2 '^x,y ^2 



^ 'TniTI*k'^Vg^x,y)\ 
(1 - iujfrnY J 



i2Tj kx^yVg(^x^y) 
(1 - iujfm.y 



(77) 



^ 1 /c^_^ 



(78) 



*h2 



-P4 = ^TTTTT I W 



1 — iuJT„ 



ir]k'^ 



1 — iuJT„ 



kxA\2yQx + kyAi2xQy 



(C + ri/3)kxky 
1 - iijfm 



[kyQx + kxQy) j , 

(79) 



iri*k'^ 



1 — iojT„ 



k^ (3C^ - 2B^D22) + 2B^ {D^k"^ - a;2j^2^ 



1 — ia;r„ 



/ {C + v/3) kxky ^ 
kxRxAi2y + kyRyAi2x z : — : [kyRx + kxRy) 

\ 1 - lUJTm 



(80) 



Q2 = ^ [k^D^ + K^)i2B^ ~ ujA^K^) ~ k^K^{2B^ + A^D^) + A^uj^K^] , = 



(81) 



Wn=Cu + By - 2D^{kyVgx + k^Vgy), 



(82) 
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{Bu , Cv 



fc2(l - iLUfm) 



1 — iuJT„ 



((1 - i2wT™)wg(j,^j^) + i2ifk^.y) , 



(83) 



«77 



1 — icjT,, 



2D^kxky - {D-y + A' 



2\ ^^mkx,y'^g[x,y) 

''' 1 - iwfm. 



7 



/ca-.y ((1 - i2u)T^)Vgi^y^^) + i27]ky^x) +i2[C+ - j k^ky 



(84) 



Qx^y — ^^x.y 



1 — iwT„ 



TmVg{y,x) [ + 2 K^kx^yVg(^x,y] 



(85) 



Wx,y =[l- 

, »(C + ^7/3) 

1 - icjTm, 



(C„ 

i — iWTm 



1 — ZWfn 



{kxBy^U ^ kyCy 2D^ ( ky x 



2i 



lOJT, 

^"^ni^x.y /7 _i_ A \ 

- ^ ^ \K.yVgx + l^xVgy) I , 



(86) 



iT-^y = ioK^ 2B^ + C22t| + + fcy522) - (3C-, + 2B^i?22) , ^22 = + {A^ + Ak" D^) Z?22, (87) 



iTmK?; 
1 — iujfrn 



(W (C22 + iK^) - kxA22 - kyB22) + (2X^^22 - C22 

2uj [uj (^22,^22) - kx,y {D^K^ + 2 {D^D22 + C22))] + kx,y (^^^22 + feyB22) (l + - ^ 
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